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Abstract. We report the detection of 22 GHz H2O maser 
emission from the red supergiant IRAS04553— 6825 in the 
Large Magellanic Cloud. It is the first known source of 
circumstellar H2 O maser emission outside the Milky Way. 
The measured flux density is comparable to that expected 
from scaling the galactic red supergiant NML Cyg. The 
peak velocity agrees with the SiO maser peak velocity. 

A near-infrared spectrum indicates 

that IRAS04553-6825 has a typical LMC metallicity. We 
argue that, possibly as a result of the low metallicity, the 
H 2 emission probably occurs near or within the dust for- 
mation radius, rather than further out as appears to be 
the case in NML Cyg and galactic OH/IR stars. 
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1. Introduction 

Maser emission from evolved stars is associated with 
dusty, oxygen rich circumstellar envelopes (CSEs), mainly 
around red supergiants (RSGs) and stars at the tip of the 
Asymptotic Giant Branch (AGB) (Elitzur 1992; Habing 
1996). The strongest maser lines originate from OH (at 
1612 MHz), H 2 (at 22 GHz), and SiO (at 43 and 86 
GHz), located progressively deeper into the CSE. Hence 
these lines are used to trace the mass-loss history and ve- 
locity structure of the CSE, through their photon fluxes 
and line profiles (Chapman & Cohen 1986; Lewis 1989, 
1990). Maser lines are also excellent tracers of the stellar 
space velocity and have been used to map the kinematics 



and mass distribution in our Galaxy (e.g. Lindqvist et al. 
1992; Sevenster et al. 1995). 

Circumstellar OH and SiO maser emission has recently 
been detected in the Large Magellanic Cloud (LMC) 
(Wood et al. 1986, 1992; van Loon et al. 1996, 1998a). 
Wood et al. (1992) argue that the outflow velocities in 
CSEs of OH/IR stars in the LMC are a factor two lower 
than in corresponding Milky Way stars, a difference they 
attribute to less efficient driving of the outflow caused by 
the lower metallicity in the LMC. However, Zijlstra et al. 
(1996) find the difference to be only 20 to 30%, from their 
analysis of the same data. Moreover, the SiO maser emis- 
sion from the LMC RSG IRAS04553-6825 (van Loon et 
al. 1996) showed that in this case, the outflow velocity 
as derived from the OH maser profile was underestimated 
by more than a factor two. The relation between outflow 
velocity and metallicity is therefore still uncertain. 

In the LMC, SiO maser emission is generally too faint 
to be detected, with IRAS04553— 6825 the only successful 
detection to date. Stellar H 2 maser emission has not been 
found outside the Galaxy before. 

Here we present the discovery of the first extra-galactic 
source of circumstellar H 2 maser emission — viz. the 
RSG IRAS04553-6825 in the LMC. We also present spec- 
troscopic results that confirm the LMC metallicity of its 
stellar photosphere. We then discuss the implications for 
the structure of the CSE by comparing with the galactic 
red supergiant NML Cyg. 
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Fig. 1. Discovery spectrum of the 22 GHz H2O maser emission 
from LMC red supergiant IRAS04553— 6825. The velocities are 
heliocentric. Also indicated are the peak velocities of the 86 
GHz SiO and 1612 MHz OH masers (arrows) and the total 
velocity extent of the 1612 MHz OH maser (bar) 

2.1. H2O radio observations 

The observations were performed on August 19, 20 of 
1997, using the 64 m radio telescope at Parkes, Aus- 
tralia. We observed the 616 — > 523 rotational transi- 
tion of ortho-H^O at a rest frequency of 22.235 GHz, 
using the 1.3 cm receiver plus autocorrclator backend. 
The 64 MHz band width and 1024 channels centred at 
~22.21 GHz yield a velocity coverage of ^860 kms -1 at 
0.84 kms -1 channel -1 . Using the Dual Circular feed we 
simultaneously obtained spectra in left and right circu- 
lar polarization. The beam FWHM is 1.4', the system 
temperature is typically 110 K, and the conversion factor 
from antenna temperature to flux density is 6.3 JyK -1 . 
The nearby sky was measured every two minutes, result- 
ing in very flat baselines that required only a very shallow 
second-order polynomial to be subtracted. 

We observed the LMC RSG IRAS04553-6825 for 
six hours on-source integration. No difference was found 
between the spectra obtained at either polarization, 
which we then averaged. The final 2 2- GHz spectrum of 
IRAS04553-6825 is presented in Fig. 1, with flux den- 
sities in mJy and heliocentric velocities in kms" 1 . The 
measured rms (1 a) is only 5.5 mJy. 

We detected a single peak of H 2 maser emission with 
a peak flux density of 68 mJy (with 20% calibration ac- 
curacy), corresponding to a 12 a detection. The peak is 
centred at a heliocentric velocity of 277 kms" 1 , and has a 
FWHM of 1.7 kms -1 . The integrated flux of the emission 
is 0.17 Jykms -1 . 

2.2. Call triplet echelle spectroscopy 

We used the 3.5 m New Technology Telescope (NTT) at 
the European Southern Observatory (ESO) at La Silla, 
Chile, on January 7, 1996, with the ESO Multi-Mode 
Instrument (EMMI), to obtain an echelle spectrum of 
IRAS04553-6825. Grating #14 and grism #4 as cross dis- 
perser were used, yielding a spectral coverage of 6000-9000 
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Fig. 2. Echelle spectra of the Ca n triplet lines in 
IRAS04553-6825 

A. The slit width and length were 2" and 4", respectively. 
The integration time was one hour. 

The data were reduced in the normal way using 
the Munich Interactive Data Analysis Software (MIDAS) 
package. The wavelength calibration was done by taking a 
ThAr lamp spectrum in conditions identical to the spec- 
trum of IRAS04553— 6825. The measured spectral resolv- 
ing power is ~ 4 x 10 4 . 

The equivalent width of the Ca II triplet lines at 8498, 
8542, and 8662 A measures the surface gravity in giants 
and supergiants, if the metallicity is known. Reversely, if 
an estimate of the surface gravity is available from knowl- 
edge of spectral type and luminosity, the Ca 11 lines can 
be used to infer the metallicity. 

The echelle spectrum of IRAS04553— 6825 around the 
Ca II triplet is presented in Fig. 2, on an arbitrary flux 
scale. The equivalent widths of the two strongest compo- 
nents are measured as W 7 8542 = 1.6 ± 0.1 A and W / 8662 = 
1.7 ± 0.1 A. Their sum is 1^8542+8662 = 3.3 ± 0.2 A. 

The Ca 11 absorption is maximum at a heliocentric ve- 
locity of ^300 kms -1 , corresponding to a redshift with 
respect to the stellar velocity of a little less than the out- 
flow velocity in the outer parts of the CSE. This is ex- 
plained by scattering of the photosphcric spectrum by an 
extended, expanding dust shell. The measured equivalent 
widths of the Ca 11 lines are not affected by the scattering 
(see Romanik & Leung 1981, and references therein). 

For IRAS04553— 6825 we adopt a bolometric luminos- 
ity of L ~ 4 x 10 5 L Q , and an effective temperature be- 
tween 2500 and 3000 K (spectral type: M7 supergiant). 
The surface gravity must then be about logg <~ cms -2 . 
If we then compare the measured equivalent width of the 
Ca II triplet with the results from Garcia- Vargas et al. 
(1998), we find that consistency is only achieved if the 
metallicity of IRAS04553— 6825 is a factor two to three 
lower than solar. This result is not very sensitive to the 
values of the adopted stellar parameters, and hence we 
conclude that IRAS04553-6825 has a typical LMC metal- 
licity (Z = 0.008). 
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3. Discussion 

3.1. Circumstellar H2O masers in the LMC 

The detection of IRAS04553-6825 shows that H 2 
masers in the LMC with photon fluxes of 10 44 s -1 are 
detected at a 5 a level within 6 hours of on-source inte- 
gration, with the current receiver at Parkes. Typical ob- 
served H2O maser photon fluxes from galactic sources are 
between 10 43 s" 1 (Lindqvist et al. 1990) and IO 44 - 45 s -1 
(Nyman et al. 1986) for OH/IR stars, IO 41 - 44 s -1 for ir- 
regular and scmircgular red variables with mass-loss rates 
- 1CT 7 M yr -1 (Szymczak & Engels 1997), IO 43 - 44 s -1 
for Mira variables (Benson & Little-Marenin 1996; Yates 
et al. 1995), and 10 44-45 s -1 for supergiants (Yates et al. 
1995). This implies that many of the mass- losing RSGs 
and AGB stars (i.e. IRAS point sources) in the LMC are 
expected to have detectable H 2 maser emission. 

3.2. Introduction to IRAS04553-6825 

With a bolometric luminosity Mbd = —9.3 mag and a 
spectral type M7, IRAS04553— 6825 was recognised to be 
the most luminous (very) red supergiant in the LMC by 
Elias et al. (1986). The progenitor mass was estimated 
M ZA MS ~ 50 Mq by Zijlstra et al. (1996). The current 
mass-loss rate of M ^few 10 -4 M Q yr _1 was estimated 
from the self-absorbed 10 fim silicate dust feature by 
Roche et al. (1993), who argue that the lower than ex- 
pected extinction in the optical indicates a disky CSE. 

IRAS04553-6825 exhibits strong OH 1612 MHz and 
weaker 1665 MHz mainline emission (Wood et al. 1986, 
1992). Comparison of the SiO maser peak velocity and 
the OH maser emission profile revealed that OH is only 
observed from the blue-shifted part of the shell. The out- 
flow velocity is 27 kms -1 , similar to Milky Way RSGs. 
(Wood et al. derived 10 kms" 1 as is suggested by the 
OH maser emission profile alone.) It is possible that the 
red-shifted emission is much weaker than the blue-shifted 
and therefore below the detection limit. Stronger blue- 
shifted emission is expected if the maser amplifies emis- 
sion from within the shell. Alternatively, an asymmetric 
and structured OH profile can also arise from bipolar out- 
flow (Chapman 1988; tc Lintel Hckkert et al. 1988). 

3.3. Similarities between IRAS04553-6825 and NML Cyg 

IRAS04553-6825 is remarkably similar to NML Cyg, a 
well known RSG in the Milky Way (Johnson 1967; Di- 
amond et al. 1984; Richards et al. 1996; Monnier et al. 
1997), in terms of luminosity, progenitor mass, spectral 
type, pulsation period, mass-loss rate, 10 feature pro- 
file, outflow velocity, CSE geometry, OH maser emission 
profile and peak flux density (scaled to a common dis- 
tance of 1 kpc), and SiO maser photon flux. The proper- 
ties of the two stars are compared in Table 1. We com- 
piled the infared photometry from the data provided in 
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Gezari et al. (1993), and use F 8flm (zeromag) = 2 x 10 -12 
Wm" 2 /im _1 = 42.7 Jy. The H 2 properties (NML Cyg: 
Yates et al. 1995) are also included in Table 1. 

The pump efficiency of SiO masers is known to be re- 
lated to the pulsation amplitude (Alcolea et al. 1990): 
Mira variables and semi-regular variables with visual 
amplitudes exceeding 2.5 mag always reach the maxi- 
mum efficiency, while variables with smaller amplitudes 
are usually less efficient. We recalculated the pump ef- 
ficiency of the SiO masers Fsio/Fs^m- The pumping in 
IRAS04553-6825 is twice as efficient as in NML Cyg, but 
a factor ~ 4 less than in Mira variables. The SiO maser 
emission of IRAS04553— 6825 has a larger ratio fint/-Fmax 
than NML Cyg (at a common resolution of 0.7 kms -1 ), 
and a few times as large as Mira variables (Alcolea et al. 
1990; van Loon et al. 1996). 

The near-infrared (JHKL) amplitudes are 0.4 mag 
for NML Cyg (Gezari et al. 1993) and 0.3 mag for 
IRAS04553-6825 (Wood et al. 1992). For NML Cyg, 
Kholopov et al. (1985) list a magnitude range of 11.19 
to 12.54 mag in R, and 17.0 to 18.0 mag in V. We have 
sparse optical photometry of IRAS04553-6825 in the pe- 
riod 1994-1997, suggesting an amplitude in V of ~ 2 mag 
(smaller in R). We conclude that the pulsation properties 
of NML Cyg and IRAS04553-6825 are similar, with pos- 
sibly IRAS04553— 6825 closer resembling Mira variables. 

From the similar properties of both objects we pre- 
dicted to detect 22 GHz H 2 maser emission from 
IRAS04553-6825 at a level of ~ 60 mJy (van Loon 
1998b). The measured peak flux density is in excellent 
agreement with this prediction. 

3.4- Location of the H 2 masers 

A striking 
difference between the H 2 masers of IRAS04553— 6825 
and NML Cyg is that in IRAS04553-6825 the emission 
peaks at the stellar restframe velocity, whereas in NML 
Cyg the line profile resembles the OH maser line profile. 

Interferometric observations of the masers around 
NML Cyg indicate that its H 2 masers are amplified ra- 
dially and that they are located in the dusty part of the 
CSE, where radiation pressure on the grains accelerates 
the outflowing matter (Richards et al. 1996). Indeed, in 
NML Cyg the strongest emission occurs at a blueshift of 
~ 20 kms -1 . 

The peak of 

H 2 maser emission from IRAS04553— 6825 falls within 
<~ 1 kms -1 of the stellar velocity (as inferred from the 
SiO velocity), and all emission is detected within a range 
of 8 kms -1 . This is similar to the situation in Mira vari- 
ables, indicating low outflow velocities in the H 2 masing 
region of the CSE (Yates et al. 1995). In contrast, NML 
Cyg is similar to the OH/IR stars. 

Cooke & Elitzur (1985) derive a scaling relation for 
estimating the inner radius r q of the H 2 masing region, 
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Table 1. Comparison between the properties of the red super- 
giants NML Cyg in the Milky Way and IRAS04553-6825 in 
the LMC. 



Bolo. Luminosity 
Distance 
Progenitor Mass 
Spectral Type 
Pulsation Period 
Mass-loss Rate 
10 /im Silicate 
(K — L) colour 
(H - K) colour 
Outflow Velocity 
CSE Geometry 
OH Profile 
OH Peak (1 kpc) 
OH Photon Flux 
SiO Profile 
SiO Peak (1 kpc) 
SiO Photon Flux 

-Fsio/ Fsfxm 

F int /F max (SiO) 
H 2 Profile 
H 2 Peak (1 kpc) 
H 2 Photon Flux 



NML Cyg 



metallicity the dust formation and the acceleration of the 
outflow occurs at larger radii because of the relative lack 
of refractory elements: in this case the difference between 
NML Cyg and IRAS04553-6825 could be a generic dif- 
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LMC. 

Acknowledgements. We would like to thank Drs. Marcus Price 
and Ian Stewart for help with the observations at Parkes. We 
Roland Gredel for help with the NTT obser- 



IRAS04553 
4 x 10 5 L Q 
50 kpc 
~50 M 
M7 

930 days a i so thank Dr. 

~5xl0 Mq yr vations at La Silla, that were performed in Director's Discre- 
emission+absorptiQfenary Time. We thank Dr. Anita Richards for helpful dis- 
1.9 mag cussions, and an anonymous referee for her/his remarks that 

1.2 mag helped improve the paper. This research was partly supported 

27 km s~ 4 by the NWO under Pionier Grant 600-78-333. O jacco agradece 

disk + bipolar a o anjinho Joana por ser para ele como ele queria ser para ela. 
blue asymm./duplicit 
1.0 kjy 
3.7 x 10 45 s" - 

peak + broad comAlcolea J.. 

270 Jy 

2.6 x 10 
j_ 

43 

15.7 km s 



,45 



stellar velocity 
170 Jy 



below which the density is higher than ~ 10 cm~ 3 and 
the maser is quenched. This was confirmed by observations 
(Yates & Cohen 1994). The collision rate is dominated 
by collisions with H 2 and does therefore not depend on 
metallicity: we therefore expect that the same relation will 
hold in the LMC. With an outflow velocity at r q of 0.8 
kms" 1 (half the FWHM of the H 2 maser peak), and 
an average particle mass of 2 x 10 -24 g, we estimate r q ~ 
4xl0 15 cm, or ~ 20 times the stellar radius. This is outside 
the region where the SiO maser occurs (Humphreys ct al. 
1996). 

We speculate that in IRAS04553-6825 the H 2 
maser originates near the dust formation radius, 
where the logarithmic velocity gradient is very high 
(d(ln(v))/d(ln(r))^l) and tangential amplification dom- 
inates over radial amplification. Part of the emission may 
also originate from inside the dust formation radius where 
outflow velocities are low. Both cases could explain the 
narrow H 2 peak coinciding with the tangcntially ampli- 
fied SiO maser peak (see also Engels et al. 1997). The H 2 
mascrs in NML Cyg are located outside of the dust forma- 
tion region, where the outflow velocities are already <T5 
kms -1 (Richards et al. 1996). 

In this model, either the maser occurs closer to the star 
in IRAS04553— 6825, or dust formation takes place further 
out (as in galactic Miras as compared to galactic OH/IR 
stars). The first possibility appears less likely if the inner 
radius is set by collisional quenching of the population in- 
version. It would be interesting to test whether at lower 
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